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@ Procddd et dtspositif penmettant ^ un modem de se synchroniser sur un transmetteur de donndes 
num^riques par vole hertzienne en presence de brouilleurs. 
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@ Le proc^d^ consists ^ estimer (9) la matrice 
de correlation ^xx. des signaux re^us du trans- 
metteur sur un ensemble de n capteurs de 
reception, d caicuier (10) la matrice de correla- 
tion inverse 1%^, d caicuier (11) des vecteurs 
d'intercorreiation fxdi entre les signaux regus 
sur I'ensemble des n capteurs et un signal de 
r^plique d(k) connu, d caicuier un crit^re de 
synchronisation (15) en effectuant les produits 
des vecteurs d'intercorreiation transposes 
Gonjugues ((i^ et de la matrice des vecteurs 
d'Interconreiation fxdi, et d comparer (18) la 
valeur du critdre obtenu ^ une valeur de seuil t) 
determinee pour placer la synchronisation sur 
I'echantllion pour lequel ia valeur des criteres 
depasse ia valeur de seuil -n. Pour estimer le 
decaiage Doppler, la sequence de synchronisa- 
tion est decoupee en tranches d'un nombre 
detennine de symboles pour effectuer des cal- 
culs de transfomnation de Fourier rapide sur 
des vecteurs d'intercorreiation detenmines sur 
cheque tranche. 
Application : Transmission radio de donnees. 
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La pr^sente invention cencerne un proc^d^ et un dispositif permettant k un modem de se synchroniser 
sur un transmetteur de donates num^riques par voie hertzienne notamment en pr^ence de brouilleurs. 

Dans les syst6mes de transmissions de donn6e8 num6riques par voie hertzienne la synchronisation est 
effectu^e. S partir de sequences de symt)oles connus du r^cepteur et Intercal^es entre les symboles d'lnfor- 
5 mation. par correlation entre le signal re9u et un signal de reference construit d partir des symtx^les connus. Le 
r6sultat est compar6 k un seuil pour d6cider de la synchronisation. Des syst6mes penmettant la synchronisation 
d'un nrtodem d partir d'un seul capteur et d'un seul r6cepteur sont par exemple connus de Tartide D.V. Sarwate et 
M.B. Pursely intitule "Cross correlation properties of pseudo random and related s^uences" paru dans IEEE Vol. 
66, pp. 593-619, May 1980. Mais les techniques mises en oeuvre qui pennettent en outre une estimation des ca- 
10 ract6ristiques du canal de transmission deviennent rapidement ineff icaces en pr^ence de brouillage. 

Lorsque le rapport signal ^ bruit n'est pas suff isant pour la prise de synchronisation par une technique 
mono-capteur, des techniques visant k accroftre le rapport signal k bruit et k permettre la synchronisation sont 
utilis^es. 

Ces techniques peuvent Stre regroup^es en deux families. Une premiere famille vise k pr^traiter le signal 
15 regu avant d'effectuer la synchronisation du signal pr^trait^. La seconde familte vise k am^liorer le rapport 
signal k bruit et k effectuer la synchronisation conjoin tement 

Dans les techniques relevant de la premiere famille il faut 6gaiement distinguer deux categories : 

- Une premiere cat^gorie met en oeuvre un pr^traitement mono-capteur de type excision de brouilleurs 
avant d'effectuer la synchronisation. Ce traitement est connu par exemple de {'article de L.B. MILSTEIN 

20 "Interference rejection Techniques in spread spectrum communications", paru dans Proc. IEEE, Vol. 76. 

pp. 657-661, June 1988. Suivant cet article on cherche k pr^dire la contribution du brouiileur, dont la 
bande est suppos^e plus petite que la bande du signal, sans pr^dire celle du signal utile. Cela permet 
d'ameiiorer les performances de la synchronisation en presence d'interf^rences bande 6troite (type 
FSK). Mais le traitement reste inefflcace en presence d'un brouillage occupant toute la bande utiiis^e 

25 par le signal utile. Cela est en partie due au fait que les techniques de pr6traitement mono-capteur ne 

permettent un traitement des signaux que dans un espace limits k I'axe des temps et d I'axe des fre- 
quences. 

- La seconde categorie concerne les techniques mettant en oeuvre un pretraltement multi-capteurs de 
rejection, precedent la prise de synchronisation. L'idee consiste k mettre en oeuvre un f iltrage spatial, 

30 discriminant le signal et les brouilleurs par leur direction d'arrivee, le temps, le spectre ou la puissance 

comme cela est decrit dans I'ouvrage de R A Monzingo, T.W. Miller, intitule "Introduction to adaptive 
arrays", 6dite par John Wiley and Sons, New-York, 1980 et la these de P. Chevalier. Intituiee "Antenne 
adaptative : d'une structure lineaire e une structure non lineaire de Volterra", Th6se de doctorat, Uni- 
versite de Paris Sud, Juin 1991 puis de correier le signal de reference avec la sortie de I'antenne. Les 
35 techniques de f iltrage spatial discriminant le signal et les brouilleurs par le temps, le spectre ou la puis- 

sance sont generatement sous-opti males, alors que celtes exploitant la connaissance de la direction 
d'arrivee du signal necessitent la mise en oeuvre d'une technique de goniometrie, coOteuse et delicate 
dans la gamme des frequences HF notamment du fait des multl-trajets. 
Pour ces raisons. les techniques de la seconde famille, visant e effectuer la rejection de brouilleurs conjoin- 
40 tement avec la prise de synchronisation, d partir de la connaissance du signal de reference, apparaissent plus 
prometteuses. Une telle technique a par exemple ete proposee par M. Compton dans son article "An adaptive 
array in a spread spectrum communication system" paru dans Proc. IEEE, vol. 66, pp. 289-298, March 1978, 
pour la synchronisation de signaux PN. k partir d'une architecture analogique integrant des "Delay Lock Loop". 
Une seconde technique f ut proposee par L.E. Brennan et I.S. Reed dans leur article intitule "An adaptive 
45 array signal processing algorithm for communications", paru dans IEEE Trans-Aero-SysL vol. AES-18, n**1, 
pp. 1 24-1 30, June 1 982 pour les transmissions numeriques et compietee par D.M. Duglos et R.A Scholtz dans 
I'artide intitule "Acquisition of spread spectrum signals by an adaptive anray", paru dans IEEE Trans-Acous- 
Speech-Signal-Proc., vol. ASSP-37, n**8. pp. 1253-1270, August 1989. Ces techniques apparaissent toutefois 
beaucoup plus couteuses que les techniques mono-capteur et sont malheureusement inapplicables telles quel- 
le les en presence des derives Doppler du signal Inhere ntes aux transmissions ionospheriques. 
Le but de I'invention est de pallier les inconvenients precites. 

A cet effet, I'invention a pour objet un precede permettant e un modem de se synchroniser sur un trans- 
metteur de signaux numeriques par voie hertzienne en presence de brouilleurs consistant k: 

- estimer la matrice de correlation Rxx de signaux regus sur un ensemble de N capteurs de reception par 
55 blocs composes d'un nombre d'echantillons d'une sequence de synchronisation transmise par le trans- 
metteur, 

- e calculer la matrice de correlation inverse R^, 
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- d calculer des vecteurs d'intercorr6lation rxd entre les signaux X(k) re9us sur rensemble des N capteurs 
et un signal de r^plique d(k) connu, 

- d calculer un entire de synchronisation effectuant le prodult scaiaire entre le vecteur d'intercorr6latlon 
txd et le vecteur obtenu par produit de I'lnverse de la matrice de correlation R;^ et du vecteur d'inter- 
corr^lation rxd 

- d comparer la valeur du crit^re obtenu h une valeur de seuil rt d6termin6e pour placer la synchronisation 
sur r^chantlllon pour lequel la valeur des entires d^passe la valeur de seuil r\, caract^risd en ce qu'il 
consiste pour estimer le d^calage Doppler inherent aux communications entre 6metteurs et r^cepteurs, 
d d^couper la sequence de synchronisation en un nombre determine de M tranches de ko symbotes, k 

effectuer sur ces tranches des calcuts de vecteurs d'intercon^eiation partielles rxdi (1=1 M) de telle 

sorte que sur chacune des tranches de ko symboles, le d^calage Doppler fasse tourner la phase d'une 
valeur inf^rieure S 90"^, d effectuer des calculs de transform^e de Fourier rapide (FFT) sur les vecteurs 
Txdi obtenus, et d calculer une valeur du crit^re de synchronisatk)n d partir des vecteurs rxd^kAf), obtenus 
en r^sultat des calculs de transform6e de Fourier rapide, et de la matrice de correlation Rxx- 

L' invention a ^gaiement pour objet un dispositif pour la mise en oeuvre du precede. 

Le precede et le dispositif selon I'invention ont pour avantages qu'lls permettent ^ la fois d'assurer la syn- 
chronisation d'un modem dans des conditions de brouillage severe et en presence de multi-trajets. d'estimer 
la derive en frequence entre ^metteur et r^cepteur et de d^tecter les instants d'arriv^e de tous les multi-trajets 
associes au signal utile, pourvu que ceux-ci soient de puissance suff isamment forte par rapport d la puissance 
de bruit de fond. 

D'autres caracteristiques et avantages de I'invention apparaitront dans la description qui suit faite au re- 
gard des figures des dessins annexes qui reprdsentent : 

- La figure 1 un schema fonctlonnel d'une antenne adaptative selon I'invention. 

- La figure 2 un exemple de mise en oeuvre d'un f litre adapts spatial d discrimination par la nrradulation 
pour la realisation d'une antenne adaptative selon I'invention. 

- La figure 3 un schema d'antenne adaptative pour la mise en oeuvre du precede de synchronisation spa- 
tiale selon I'invention. 

- La figure 4 un graphe pour illustrer la methode d'estimation de la matrice de correlation de signaux pro- 
venant de N capteurs. 

- La figure 5 un graphe pour illustrer le decoupage des symboles de la sequence de synchronisation par 
tranche pour le calcul de vecteur d'intercorreiation en presence de decalage Doppler. 

- figure Bun graphe montrant revolution de critere de synchronisatton mis en oeuvre par I'invention 
en presence d'un decatage Doppler. 

- La figure 7 un organigramme retra^ant ['ensemble des etapes necessaires d la mise en oeuvre du pro- 
cede de synchronisation selon I'invention. 

Le precede de synchronisation mis en oeuvre par I'invention fait appel aux techniques connues de syn- 
chronisation mono-capteur pour les appliquer d la realisation d'antennes adaptatives sur replique. 

Selon la technique de synchronisation mono-capteur, le critere de synchronisation C(n) est obtenu en cal- 
culant la correlation y(n) entre K symboles connus dans une sequence de synchronisatton et les K echantillons 
qui suivent I'arrivee d'un echantillon n sur le capteur, en normalisant les resultats par la puissance p(n) du signal 
repu. Ces calculs sent effectues en appliquant les relations : 



y(n) = 7i;x(k+n)d*(k) (1) 
k=i 



1 ^ 

P(n) = 7E'«(k-»'n)x«(k+n) (2) 

C<n)=l^ (3) 

oil: 

- x(k) est le signal issu du capteur utilise, 

- d(k) est le signal de replique connu (K symboles places en debut de trame), 
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- * repr6sente reparation de cxDnjugaison sur les nombres complexes. 

Ces calculs sent r6it6r6s pour cheque 6chantillon n, et lorsque le critfere C(n) d6passe un certain seuil r\, 
dependant de la probabilit6 de fausse alarme, la synchronisation est d6cld6e. L'Instant de synchronisation est 
d6termin6 par le rang d'6chantillon n© pour lequel C(no)> r\. La normalisation par la puissance p(n) permet de 
choisir un seuil independent du niveau de signal regu. 

Dans ce proc^dd la synchronisation est obtenue lorsque le rapport signal d bruit est suff isant Ses per- 
fonmances dependent en partie du nombre d'6chantillons pris en compte pour la correlation. Cependant en 
presence de brouilleurs forts, la synchronisation ne peut dtre effectude. 

L'id^e mise en oeuvre dans I'invention conslste alors d effectuer la synchronisation, en sortie d'un traite- 
ment d'antenne adaptatlve, permettant la rejection des brouilleurs. 

Le schema fonctionnel d'une antenna adaptative est represents d la figure 1 . celle-ci comporte un rSseau 
1 de N capteurs 1 1 ... 1 n, un rSseau de f litres adaptes 2 (constitues par des circuits de ponderatlons amplitude- 
phase, des f litres d r6ponse impulsionnelle f inie RIF ou inf inie Rll) couple ^ un processeur de commande adap- 
tatif 3. Le processeur 3 calcule les coefficients Wi ^ du reseau de f litres 2 pour I'adapter au signal utile 
et aux brouilleurs, en Tabsence de toute connaissance a priori sur les brouilleurs et k partir d'une Infonmation 
minimale sur le signal utile. Les coefficients d Wn pemnettent de former un trou dans le diagramme de 
rayonnement du reseau d'antenne dans la direction des brouilleurs, tout en maintenant un certain gain dans 
la direction du signal utile. 

Differents types d'information sur le signal utile sont exploites par Tantenne adaptative. Ceux-ci sont d'or- 
dre spatial lorsque la direction d'anrivee du signal utile est connue, d'ordre temporelle lorsque les Instants de 
presence et d'absence de signal utile sont connus par example en evasion de frequence, relatifs k la forme 
d'onde lorsqu'une information sur la modulation du signal utile est disponible pour disposer d'une replique du 
signal utile en intercalant par exemple des sequences connues dans la trame utile, d'ordre spectral lorsque 
le signal utile occupe une certaine bande de frequence et que les brouilleurs occupent une bande plus impor- 
tante, ou statistiques pour separer les dif ferentes sources arrivant sur le reseau en exploitant par exemple les 
statistiques d'ordre 4 des signaux en supposant que les differents signaux sont statistiquement independents. 

Dans le cas de I'antenne adaptative representee e la figure 2, 1'information sur le signal utile est constituee 
par les sequences de symbojes connus intercaiees dans la trame entre les sequences de symboles d'infor- 
mation. Les f litres places en sortie de chacun des capteurs 1 1 ^ 1^ sont formes respectivement par des circuits 
de ponderation amplitude-phase Wi d Wn des signauxxi(t) e Xn(t) fournis paries capteurs. Les signaux obtenus 
en sortie des circuits de ponderation sont appliques d un circuit de sommation 4. 

Les caracteristiques du reseau de f litres 2 sont calcuiees en minimlsant une Erreur Quadratique Moyenne 
(EQM) estimee I entre un signal de replique d(k) et le signal y(k) sortant de I'antenne de la faffon schematises 
e la figure 2 ou les elements homologues^ ceux de la figure 1 sont reper6s avec les mfimes references, \ est 
definle par la relation 

K K 

i = Z|y(k) - cl(k)p = 2 |w*X(k) - d(k)f (4) 
k=1 k=l' ' 

oCi: 

- y(k) est la sortie de I'antenne adaptative et est definle par 

y(k) = WX(k) (6) 

- W est le vecteur poids def inissant I'antenne adaptative, W etant le vecteur transpose conjugue corres- 
pondanL (le signe t en exposant def init t'operateur correspondent.) 

- X(k) est le vecteur forme par les signaux re^us sur les capteurs. 

- d(k) est le signal de replique connu. 

La difference y(k) - d(k) est obtenue sur la figure 2 e la sortie d'un circuit soustracteur 5. 
Le vecteur poids W est calcuie par la formula connue de Wiener : 

W = R^^fxd (6) 

avec : 



Rxx = ^ix{k).xHk) (7) 
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matrice de correlation des signaux provenant des capteurs 1 1 1 n- (X^k) 6tant le vecteur transpose conjugu^ 
du vecteurX(k)). 



rxd = ;^i:x(k)d*(k) 



vecteur d'intercorr^lation entre X(t) et d(t) (8). 

10 L'antenne adaptative ainsi utilis^e conduit d una maximisation du rapport slgnal/(brult + interference) 
connu sous la designation SNIR abr6viatlon anglo-saxonne de "Signal to Noise + Interference Ratio" obtenu 
en sortie de l'antenne, avec une possibility de rejection pouvant aller jusqu'd N-1 brouilleurs, et est appeiee 
Filtre Adapte Spatial (FAS). 

recherche de la synchronisation est effectuee selon Tinvention en sortie de l'antenne adaptative en 

15 discriminant le signal selon sa modulation, et en utilisant le critere de synchronisation mono-capteur. Le sche- 
ma de principe de la synchronisation spatiale est represente d la figure 3 ou les elements homologues d ceux 
des figures 1 et 2 sont representes avec les memes references. Le critere de synchronisatton devient, en uti- 
lisant les notatk)ns precedentes : 

20 

C(n) = l|;y(k + n)d(k) (9) 

25 y(k+nji represente la sortie de l'antenne adaptative conduisant d une minimisation de I'erreur quadratique 
moyenne |(n) estim^e sur K §chantillons entre la r^pllque d(k) et le signal arrivant pris d partir de I'tehantillon 
n, suivant la relation 



K 

^(n)= i;|y(k+n)-d(k)f avec y(k+n) = W(n)*.X(k+n) (10) 
k=1 



avec : 

35 W(n) = Rxx-Hn).fxd(n) (11) 

Les diff^rentes correlations sont calcul^es sur les K echantillons qui suivent I'arrivee d'un echantillon n 
suivant les relations : 



1 

.Rxx{n) = -5:x(k + n)xHk+n) (12) 
K 

rxd(n)-^2;x(k + n)d*(k) (13) 



Lorsque rechantillon n correspond e la position de synchronisation, les poids W(n) definissent le filtre 
50 adapte, et y(k-i-n) contient alors le signal utile debarrasse de la contribution des brouilleurs. Dans le cas contrai- 
re. aucun signal n'est correie avec la replique d(k) et y(k+n) ne contient plus alors que la contribution du bruit 
de fond. 

L'expression du critere de synchronisation spatiale se simpllf ie en ecrlvant : 



C(n) = 1 2 W(n)*X(k + n) d*(k) = W(n)* fl ^X{k + n) d*(k) (14) 

k=i k=i 
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soit 

C(n) = ^Jn)Rxx(nKxd(n) (15) 
De m6me que pour le cas mono-capteur, ce calcul est r6it6r6 pour chaque 6chantillon n. et lorsque le crit6re 
C(n) d^passe un certain seuil ^, dfependant de ta probability de fausse alanme. la synchronisation est d6cid6e 
(ia normalisation par la puissance est ici obtenue par rinterm6diaire de R,^^), 

L'instant de synchronisation est alors d6termin6 par te rang de l'6chantillon pour lequel C(no) est sup6rieur 
d une valeur de seuil d4termln§e r\. 

Les calculs correspondant sont effectu6s sur la figure 3 par le bloc 7 de calcul du entire C(n) et par le 
d^tecteur de seuil 8, ^ partir de la sequence de synchronisation d(t) et du signal y(t) corr6l6s par le corr6lateur 
6. 

Cela n6cessite pour chaque 6chantillon, d'une part de calculer N correlations entre le signal de r6plique 
d(k) et les signaux x(k) re^us par les N capteurs (soit N fois les calculs effectu6s pour la synchronisation mono- 
capteur), et d'autre part d'estimer puis d'inverser ta matrice d*intercorr6Iation des signaux capteurs. 

Au cours de ces calculs la synchronisation et le calcul du filtre adapts 2 sont effectu6s oonjointement. 
Lorsque la valeur du crit^re d6passe le seuil ti, calcul6, la synchronisation est d6cid6e et te vecteur poids d6- 
f inissant le fittre 2 s'en d6duit facilement par la fonmule connue de Wiener. Le vecteur poids W peut ensuite 
§tre utilise pour f iltrer les 6chantlllons contenant rinformation. 

En supposant que la r^plique d(t) est de puissance 1 , les correlations tendent vers les valeurs suivantes, 
lorsque le nombre d*echantitlons K devient tr^s grand : 

Rxx = R + KaS S< 

ou: 

- R est ta matrice de correlation des Interferences et du bruit de fond, 
2s - S est la matrice du vecteur directeur du signal utile, 

- 7Cs est la puissance du signal utile regu sur un capteur. 

Le critere de synchronisation qui converge alors vers une valeur C© s'ecrit encore: 

- ^ S _ SNIRopt 
° 1 + 7c,StR-iS " 1 + SNIRopt ^ ^ 
30 OU SNIRopt est le SNIR asymptotique obtenu en sortie du filtre spatial adapte 2. 

Lorsque le SNIR est fort, le critere tend vers 1, lorsqu'il est falble, it tend vers 0. L'importance du choix du 
seuil Ti, fixe la probabitite de detection de la synchronisation. Le seuil depend naturellement de la probabilite 
de fausse atanme souhaitee. 

Le calcul du critere de synchronisation C, de fxd et de Rxx est effectue pour chaque position de la syn- 
35 chronisation et notamment la matrice Rxx est estim6e et inversee egalement pour chaque position de la syn- 
chronisation, la puissance de calcul demandee est eievee. Mais it est possible de reduire la puissance de calcul 
en estimant la matrice de correlation Rxx moins souvent sur un intervalle de temps plus important (K' echan- 
tillons au lieu de K avec K>K de la fagon representee e la figure 4) et en utilisant cette estimation de Rxx sur 
toutes les positions de synchronisatijpn incluses dans cet Intervalle de temps. 
40 Dans ces conditions la matrice Rxx s'ecrit: 



20 



+ 8 (17) 



50 oil: 

- B(tk) est la contribution du bruit et des brouilleurs sur X(tK), 

- S(tk) est la contribution du signal utile sur X(tO, 

- Ko est rinstant d'apparitlon du signal utile sur le palier, 

- e est un terme prenant en compte rintercorreiation entre le signal B(t) et le signal utile S(t), qui tend vers 
55 0 lorsque K devient grand. 

Pour la position de synchronisation testee : 



45 



Rxx = 



ZB(tk)B*(tK)+ is(gsHtK) 



k=1 



k=Kn 
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Ko+K 
Z X(tk)d(k) 
k=Ko 



(18) 



10 
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Dans les conditions asymptoliques (avec K»1), lorsque le signal utile est synchronis6 avec la r6plique 
d(Q> ^xx 6t Txd se simpiif lent at s'^crivent : 

Rxx^R"- \ SS^ (19) 



A 

rxd^ 



,^^/^^S (20) 

Le crit6re obtenu en estimant la matrice de correlation Rxx sur K' 6chantlllons poss^de comme valeur 
asymptotique, apr^s un calcul utilisant le lemme d'inverslon matricielle: 

^ SNIRopt 



1 + Kjl^jLl SNIRopt 



(21) 



K' — Kn + 1 

Comme ^ < 1 , alors Co > Co 



La valeur du entire apparalt dans ces conditions augment6e d'un facteur dependant de Tinstant d'appa- 
rition du signal utile sur la frequence consid6r6e, ce qui entraine une amelioration de la prise de synchronisa- 
tion, au prix d'une puissance de calcul r^duite. 

D'autre part, le nombre K* ne doit pas 6tre choisi trop grand, sa valeur est d6termin6e pour que sur 
K' echantillons, Tenvironnement de brouillage puisse Stre consider^ comme stationnaire. 

Le seuil ^ prendre en consideration depend de la probabilite de fausse alarme (Pfa). La (Pfa) est la pro- 
babilite qu'un bruit soit detecte comme une synchronisation. Le calcul de la valeur du seuil a lieu en calcutant 
Pfa=P(C>Ti) lorsque I'antenne regoit uniquement du bruit. 

En presence de bruit le critere de synchronisation s'6crit: 



30 



C = 



K k=i j\r k^i J ; 



(22) 



^ En supposant que est un bruit gaussien complexe de puissance a^, independent entre les capteurs, 
lorsque K est suff isamment grand, le terme du milieu de la relation (22) devient : 



40 



17 Z ^k Br 
k=1 



£a+2id (23) 



45 



(Id etant la matrice identity ou unitaire) 

le critere de synchronisation peut alors fitre approxim^ de la fa^n suivante : 



C = 



ZBkC^jf|;BkC^*l (24) 



50 



Soit: 



55 



1 



N 



1^2 „2 ^ 

^ i=1 



K 

ZBk 
k=1 



2^ 



K o i=i 
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avec: 



K 

ZBl,C^* = R| + jIi (26) 
k=1 

Bk repr^sente le bruit obtenu sur ie capteur i, pour I'^chantillon k. 

Ri + jl, peut 6tre approxim6 par un bruit gausslen complexe de variance K o^. Dans ces conditions R| et l| 
sont chacun des bruits gaussiens r^els de variance 



(ko^) 



En notant T| = Rj ou l|, C devient : 



c=i;rTE(T|^) (27) 

^ 1=1 



K 



ou Ti est une variable al^atoire qui suit une loi nomnale centr6e de variance 1-a variable al6atoire Z| ' 
V2 

Ti suit done une lol nonmale centr6e de variance 1. 



VKa2 

C devient : 



1 2N 

C = 2^E2f (28) 



Comme la suite des Z) peut §tre consid^r^e comme une suite de 2N variables al^atolres Ind^pendantes 
gausslennes de nnoyennes nulles et de variances 1, la variable al6atolre 

2N 

1=1 

suit une lol du ^ 2N degr^s de liberty. 
La probability d'avolr C ^ n s'^crit ators : 



p(c^n) = p 



f2N 

u=i 



= Pfa (29) 



Le tableau sulvant donne diff^rentes valeurs de seulls r\ obtenues par recherche dans une table donnant 
P(2^) pour une loi du ^ 2N degr6s de liberty et pour Pfa = 10^. 
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Tl = 0,26 

T, = 0,23 
T] = 0,21 
Ti = 0,18 

Comme entre r6misston et la r6ception. le signal transmis peut 6tre effects d'un certain d^calage en fre- 
quence, ou d^calage Doppler. du essentiellement aux derives en frequences des oscillateurs locaux des 
synthetiseurs d remission et h la reception, ainsi qu'aux conditions de propagation ionospherique, il est ne- 
cessaire de compenser ce d6calage. Le risque est s'il n'est pas compense. d'empficher toute prise de syn- 
chronisation. 

En appelant Afo, le d6calage Doppler entre le signal repu par I'antenne et le signal attendu et en consid6rant 
la position de synchronisation, le vecteur signal regu peut s'ecrire avec les notations des paragraphes sul- 
vants : 

X(t) = V7^^d(t)e?l»-AVts + B(t) 
La matrice de correlation estimee peut alors s'ecrlre : 



N = 5 


=> 


A = 42 


N = 4 




A = 37 


N = 3 




A = 33 


N = 2 




A = 28 



K 

RXX = Z [V5^ d(t^)e2.j.n.Afo.t, g + B(tk)][^ ^^^2:i.nAfo\ ^ + B(1k)T (30) 
k— 1 J 



Rxx = e + 1 S[«sid(l|c)f S S» + B(tK) B(tk)»] (31) 

ou 8 est une matrice prenant en compte les Intercorreiations entre la contribution du signal utile et le bruit b(t). 
Lorsque K»1, dans les conditions asymptotiques, Rxx peut s'approximer par : 

Rxx = «»SS+ + c2|d (32) 

La relation (32) montre que le decaiage Doppler ne modif ie pas le calcul de Rxx- En effet, le decalage 
s'eiimine naturellement lorsque le signal utile est multipiie parse quantite conjuguee. 

Le vecteur d'intercorreiation peut s'ecrire : 

- K - 

^Xd = Z[V^ d(tk)e^i-^^'^ S + B(t,)] d-(tk) (33) 



^^Xd = B" + i^2:[v^|d(tK)fe2J-^fc-tKs] (34) 

ou 8' est un terme prenant en compte rinterconreiation entre le bruit B(t) et la replique. Lorsque K»1, rxd peut 
s'approximerpar: 




S (35) 



La relation (34) montre que le decalage Doppler modif ie I'estimalion du vecteur d'intercorreiation rxd et 
que du fait que la phase du tenme en exponentiel compiexe peut tourner de plusleurs fols 360° sur rintervaile 
de temps utilise pour le calcul de la correlation, I'integration ne peut plus avoir lieu de fagon coherente ce qui 
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10 



IS 



emp6che le crit^re de synchronisation multlcapteurs de converger. Le processus suivant permet d'estimer le 
d^calage Doppler et de le connpenser pour rendre possible la prise de synchronisation. 

Selon ce processus, les K symboles de la sequence de synchronisation sont d^coup^s de la fa^on repr^ 
sent^e d la figure 5 en M tranches de Kq symboles (M Ko = K), de telle sorte que sur chacune des tranches 
de Ko symboles, le d^calage Doppler fasse tourner la phase du signal regu par Tantenne d'une valeur inf6rleure 
i 90°, af in de pemriettre une integration coh6rente de rxd sur chacune des tranches. 

Si rxdt d6signe la valeur de l'intercorr6lation estim6e sur la tranche i, celle-ci s'6crit : 



^Xdi = 



^ k=iKo+1 



S i = 0....,M-1 (36) 



En mettant en facteur dans la relation (36) un terme dependant uniquement de Ko, la longueur de la tranche 
6l6mentaire consid6r6e. et du d6calage Doppler. la relation (36) devlent : 



20 



'Xdi = 



i^|>/5irK(k4iKo)T)fe'^-f' 



2M.-^JKo 



•S i=0 M-1 (37) 



25 



30 



Ce qui peut encore s'^crire : 



rxdi B X, e?l»-^iK, .s i = 0 M - 1 (38) 



F, 



La relation (38) montre que les vecteurs rxd sont des vecteurs, 6chantillonnds k -r. Leur transform^e de 

Ko 

Fourier discrete est alors : 



35 



M-1 



rxd(Af)s2:0<die (39) 
1=0 



40 



M-1 2i7u^^^iK» 
i=0 



(40) 



La relation (40) montre que lorsque Af = Af©, le d^calage Doppler est compens6, ce qui rend maximum la 
valeur du entire de synchronisation rxd(AO^xx*^rxd(Af). 

Afin d'estimer Af, c*est-^-dire afin d'estimer le d6calage Doppler qui conduit d la valeur maxi male du crl- 
^ t6re, et afin de r6duire le coOt de calcul correspondant. une transform6e de Fourier rapide est effectu6e sur 
les vecteurs rxdi (le nombre de tranches ^l^mentaires M ^tant choisi pourformerune puissance de 2) : ce calcul 
donne: 



50 



M-1 -2in.p^.ik M-1 

i=0 



k = 0.....M-1 (41) 



oCi Afi est la resolution Doppler telle que : 

M 

Si la resolution est insuff isante. une amelioration d'un facteur 2 peut Stre obtenue en calculant une trans- 
formde de Fourier rapide, sur les vecteurs rxoi "decaiees** en frequence de la fa^on survante : 
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k = 0 M-1 (43) 



L'op^ration consiste d multiplier chacun des vecteurs d'intercorr^lation txdi par le terme er o»vm) et d effec- 
tuer ensulte un calcul de transform^e rapide sur les vecteurs atnsi obtenus. 

Le crit^re de synchronisation est ensuite estim^ sur chacune des cases de calcul des transform^es de 
Fourier prdc^dentes par la relation : 



c(Kf) = fi.[Kf]R;J^kf] k = 0....,2M-1 (44) 

Bien ^videmment, si la resolution est toujours insufflsante. ii est possible de i'am^iiorer en utilisant le 
m§me processus pour caiculer des vecteurs 

oCi n est un entier queiconque. 

Un test sur la case de la transform^e correspondant au maximum du crit^re de synchronisation est effeo 

tu6. Si le seuil est d^pass^ pour une valeur ko ^ donn6e du d^calage Doppler, 11 est alors possible d'aff Iner 

le calcul du d^calage Doppler (Indispensable pour la demodulation) en calculant la valeur du critSre pour dif- 

ferentes valeurs de d^calage Doppler entourant ko En effet. d Tissue de cette etape, la resolution sur Tes- 

timation du d^calage Doppler est Pour Tameiiorer, il faut proc^der par dichotomie : d chaque iteration et 

de la fa^on representee ^ la figure 6, la valeur du critere pour la valeur de decalage Doppler situee au milieu 
des deux vateurs de decalage Doppler donnant les valeurs maxima du critere lors de ['iteration precedente 
est calcuiee. La recherche par dichotomie est arretee des que la resolution souhaitee sur i'estimation du de- 
calage Doppler est obtenue. Apres n iterations de Talgorithme de recherche par dichotomie, ta resolution, qui 

etait initialenf>ent de^, est : 

^ = -^ (45) 

LenchaTnement des etapes de calcul precedente pour la mise en oeuvre du procede selon invention d 
Taide de moyens de traitement appropries formes par exemple par des processeurs de traitement du signal 
ou tout dispositif equivalent connu est figure sur Torganigramme de la figure 7. 

II consiste suivant I'etape 9 e estimer la matrice de correlation Rxx suivant la relation (7). Cette matrice 
est estimee sur un bloc contenant un nombre d'echantillons K' superieur au nombre d'echantillons de la se- 
quence de synchronisation. 

L'estimation de Rxx est utilisee ensuite pour toutes les positions de synchronisation Induses dans le bloc. 
Puis ta matrice de correlation inverse est calcuiee e retape 1 0. rxdi qui est def ini par la relation (8) est cal- 
cuie e retape 11. Pour chaque echantilton conrespondant d une position de synchronisation P dans la suite 
des echantillons P = 1 e P = IC-K, rxut est calcuie pour i = 1 e M sur M blocs de Kq symboles comme def ini d 
la figure 5 suivant la relation 




fxdi = Tr-2X(P+K + iKo) d (K + iKo) (46) 
•^0 K=1 
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Des calculs de transform6e de Fourier Rapide sont effectu6s ^ I'^tape suivante 12 sur un nombre d6ter- 
min6 M (M = 8 par exemple) de vecteurs fxdi suivant la relation (41). 

Aux 6tapes 1 3 et 14 un calcul de vecteur Rxdi d6cal6 et un calcul de transfbrm6e de Fourier rapide suivant 
la relation (43) sont ex^cut^s pour am^tiorer ta resolution par un facteur de 2. 

Le crit^re de synchronisation est calculi k l'6tape 15 en application de la relation (44). 

La comparaison du crit^re ou seuil r\ est effectu^e k I'^tape 16. Si d T^tape 16 la valeur du critdre 




est sup6rieure ou ^gale au seuil. la position de synchronisation est d6finie d r6tape 17 sur la position P de 
l*6chantillon correspondent L'estimation du d6calage Doppler peut ensuite 6tre aff in6 k i'6tape 1 8 en calculant 
h nouveau la valeur du crit^re de synchronisation pour des d^calages Doppter entourant le d^calage 

M) 

pr6c6dent. Cette 6tape peut §tre suivie par une 6tape 19 d*estimation de la position de synchronisation pour 
les autres trajets. 

La recherche de synchronisation pr6c6demment d6crite est effectu6e pour chacun des 6chantillons de 
signal num^ris^. En pratique, pour am^liorer la precision de la position de la synchronisation il est preferable 
d'effectuer un sur6chantiilonnage du signal d'un facteur determine L, L etant egal k 4 par exemple. Les trai- 
tements de I'organigramme de la figure 7 et d6crits ci-dessus sont encore utilisables sur le signal ainsi sur6- 
chantiilonne, mais il est possible d'utiliser une variante dans te precede, af in de reduire la puissance de calcul 
necessaire au traitement Dans ces conditions, aprfes avoir calcuie pour cheque echantillon le vecteur fxdi pour 
i = 1, ... M sur M biocs de Kq symboies suivant la relation (46) un crit6re de synchronisation Q est calcuie sur 
chacun des M blocs par la relation 

= i^^F^fxdi (47) 
Puis un critere somme de Tensemble des Q est calcuie suh^ant la relation : 



Ci(P)-2Ci (48) 
i=1 

Ce critere correspond e une integration non coherente des Ci (contrairement e ce qui est fait pour le calcul 
par FFT ou on integre de fagon coherente), ce qui se traduit par une remontee du crit6re sur les echantillons 
de bruit 

Le calcul du critere a lieu par blocs de L echantillons consecutifs. Sur chacun de ces blocs, rechantillon 
qui donne la valeur maximum de Ci(P) est seiectionne pour effectuer le calcul par transform6e de Fourier ra- 
pide necessaire k l'estimation du decalage Doppler comme dejd decrit aux etapes lid 19 de la figure 7. 

Revendlcatlons 

1 . Precede permettant un modem de se synchroniser sur un transmetteur de signaux numeriques par voie 
hertzienne en presence de brouiileurs consistent d : 

- estimer (9) la matrice de correlation Rxx de signaux regus sur un ensemble de N capteurs de recep- 
tion (1i ... 1n) par blocs composes d'un nombre d'echantillons d'une sequence de synchronisation 
transmise par le transmetteur, 

- d calculer (1 0) la matrice de con-eiation inverse R^^, 

- k calculer (11) des vecteurs d'intercon-eiation fxd entre les signaux X(k) regus sur I'ensemble des N 
capteurs et un signal de replique d(k) connu, 

- e calculer un critere de synchronisation (7) effectuant te produit scalaire entre le vecteur d'intercorre- 
lation rxd et le vecteur obtenu par produit de Tinverse de la matrice de correlation R;^ et du vecteur 
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d'intercorr^lation rxd 

- d comparer (8) la valeur du critfere obtenu ^ une valeur de seuil ti d6tennin6e pour placer la synchro- 
nisation sur r^chantlllon pour lequel la valeur du crit^re d^passe la valeur de seuil ti, caract^ris^ en 
ce qu'il consiste pour estimer le d6calage Doppler inh6rent aux communications entre 6metteurs et 
r6cepteurs, ^ d6couper la sequence de synchronisation en un nombre d6termin6 de M tranches de 
ko symboles, k effectuer (11,13) sur ces tranches des calculs de vecteurs d'intercorr6lation partielles 
fxdi (i = 1. .... M) de telle sorte que sur chacune des tranches de ko symboles, le d6calage Doppler 
fasse tourner la phase d'une valeur inf^rieure d 90°, ^ effectuer (12, 14) des caiculs de transfomn^e 
de Fourier rapide (FFT) sur les vecteurs fxdi obtenus, et ^ calculer une valeur du crit6re de synchro- 
nisation d partir des vecteurs rxd(kAf)j^ obtenus en r6sultat des calculs de transform6e de Fourier ra- 
pide, et de la matrice de correlation Rxx- 

Proc6d6 selon la revendication 1 caract6ris6 en ce que Testimation de la matrice de con-elation Rxx est 
obtenu par la relation 



ou p est la position du premier echantlllon S partir duquet la correlation est effectu^e entre le vecteur X 
et la replique d. 

Precede selon la revendication 2 caracterise en ce qu'il consiste ^ choisir la case de calcul de la trans- 
formee rapide qui donne une valeur maximum du critere de synchronisation et d comparer la valeur maxi- 
mum obtenue ^ une valeur de seuil determinee pour decider de la synchronisation. 

Precede selon Tune quelconque des revendications 1^3 caracterise en ce qu'il consiste pour obtenir une 
meilleure resolution sur le decalage Doppler et une valeur plus importante du maximum du critere, d cal- 
culer (1 3) des vecteurs rxdi decaies en frequence et d effectuer un calcul de FFT sur les vecteurs Txdi de- 
caies de fa^on d obtenir la vajeur du critere avec une meilleure resolution sur les frequences Doppler. 

Procede selon Tune quelconque des revendications 1^4 caracterise en ce qu'il consiste e surechantil- 
lonner les signaux regus.sur les N capteurs par rapport au rythme symbole, pour ameiiorer la resolution 
sur la position de synchronisation. 

Precede selon I'une quelconque des revendications 1 e 5 caracterisee en ce qu'il consiste e effectuer le 
calcul de la matrice de correlation Rxx sur un nombre d'echantillons K* superieur au nombre K d'echan- 
tillons de la sequence de synchronisation afin de diminuer la puissance de calcul. 

Precede selon I'une quelconque des revendications 1 d 6 caracterise en ce qu'il consiste pour diminuer 
la puissance de calcul et apres avoir sur6chantillonne d'un facteur L par rapport au rythme symbole, e 
effectuer une etape de pretraitement consistent e calculer pour cheque echantillon p (p = 1 e L) composant 
un symbole les valeurs C| du critere de synchronisation sur chacune des M tranches de ko symboles et 
d calculer la valeur du critere en effectuant la somme des criteres C| obtenus pour retentr dans la suite 
du traitement I'echantillon donnant la valeur maximale du critere C = £C|. 

Dispositif pour la mise en oeuvre du precede selon I'une quelconque des revendications 1 e 7 du type 
comportant, un reseau de N capteurs (1i...1n) de signaux de reception couple ^ un reseau defiltres adap- 
tes (2) comportant des circuits de ponderations amplitude phase commandes par un processeur de 
commande (3), le processeur de commande (3) etant programme pour 



Rxx=4Zx(k).X(;) 



et le vecteur d'intercorreiation est estime suivant la relatbn 




13 



EP 0 665 665 A1 



estimer la matrice de correlation Inverse R;^ des signaux regu X^k) et calculer les vecteurs d'inter- 
correiatlon rxd entre les signaux X^k) regus etie signal de r6pliqued(K)Caract6ris6 en ce qu'il comprend 
des moyens de traitement (5...19) pour calcuier un crit^re de synchronisation (7) en effectuant un 
produit scalaire entre le vecteur d'lntercorr^latjon Txd et le vecteur obtenu par le produltde I'inverse 
de matrice de correlation et du vecteur d'intercorr§Iation fxd» 

comparer (8) la valeur du entire obtenu ^ la valeur de seuil r\ pour placer la synchronisation sur 
i'^chantlllon pour lequel ta vaieur du crlt^re d^passe la valeur du seuil r\, 

et estimer le d^calage Doppler inherent aux communications entre ^metteurs et r^cepteurs, en d^- 
coupant la sequence de synchronisation en un nombre determine de M tranches de ko symbotes, 
en effectuant (11, 13) sur ces tranches des calcuis de vecteurs d'intercorretatlon partielles fxd^ (i - 

1 M) de teiie sorte que sur chacune des tranches de Ko symboles le d^caiage Doppler fasse tour- 

ner d la phase d'une valeur inferieure d 90'', en effectuant (12, 14) des calcuis de transformee de 
Fourier rapide (FFT) sur les vecteurs Pxd, obtenus, et en calculant une valeur du crit^re de synchro- 
nisation e partir des vecteurs rxd(KAf^, obtenus en r^sultat des calcuis de transform^e de Fourier ra- 
pide et de la matrice de correlation Rxx- 
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Calcul du critere 
C(n)=^^ y(k+n)d*(k) 

C(n)=r*j(n)R;]j(n)rxd(n) 



Synchronisation si : 

C(n) >r| 
(ri dependanf de la Probabilife 
de fausse alarme desiree) 



Processeur 



W(n)=R;;;(n)rxd{n) 




FIG.3 
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Estimation de R^^^ 
sur K' echantillons 
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-10 



P=1 



Calcul de r^dj 
i=1 ,M= 2"" 
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